Condition for emergence of complex eigenvalues 
in the Bogoliubov-de Gennes equations 
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The condition for the appearance of dynamical instability of the Bose-condensed system, charac- 
terized by the emergence of complex eigenvalues in the Bogoliubov-de Gennes equations, is studied 
analytically. It is concluded that the degeneracy between a positive-norm eigenmode and a negative- 
norm one is essential for the emergence of complex modes. Based on the conclusion, we justify the 
two-mode approximation applied in our previous work [E. Fukuyama et at, Phys. Rev. A 76, 043608 
(2007)], in which we analytically studied the condition for the existence of complex modes when the 
condensate has a highly quantized vortex. 
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I. INTRODUCTION 

Experiments of the Bose-Einstein condensatesfBECs) 
of neutral atomic gases, first realized in 1995 [J 0, [1], 
have been offering chaUenging subjects in theoretical 
foundations of quantum many-body problem. Among 
others, some unstable phenomena of the BECs are very 
interesting, since to formulate unstable quantum many- 
body systems is an open problem. Observed examples 
of such phenomena are the split of a doubly quantized 
vortex into two singly quantized vortices jj^], the decay 
of a condensate flowing in an optical lattice [5| , and the 
decay of the initial configuration in a quenched ferromag- 
netic spinor BEC [6|]- These phenomena are observed at 
very low temperatures where the dissipative mechanism 
brought by the thermal cloud is negligibly small. Such 
instability is called "dynamical instability," and is dis- 
tinguished from the "Landau instability," in which the 
thermal cloud plays a dissipative role [7|. 

For the theoretical investigations of the instability, the 
Bogoliubov-de Gennes (BdG) equations are em- 

ployed. The BdG equations follow from linearization of 
the time-dependent Gross-Pitaevskii (TDGP) equation 
[ll[ and determine the excitation spectrum of the con- 
densate. The BdG equations have complex eigenvalues 
for some cases (see below) , and the presence of the com- 
plex eigenvalues is interpreted as the sign of dynamical 
instability. This instability is associated with the decay 
of the initial configuration of the condensate and can oc- 
cur even at zero temperature. On the other hand, the 
Landau instability, which is characterized by the nega- 



tive eigenvalues and in which the thermal cloud drives 
the system towards a lower energy state in a dissipative 
way, is impossible at very low temperature. It is reported 
that the BdG equations have complex eigenvalues in the 
cases where the condensate has a highly quantized vortex 
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12 , 13 1 , where the condensate flows in an optical lattice 



15| , and where the condensate has gap solitons 16 1 , 

18|. 
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and in the case of the multi-component BECs 
It is suggested that some degeneracy between a positive 
mode and a negative one in the BdG equations is nec- 
essary for the emergence of complex eigenvalues [ij, [l3| ■ 
The analytical investigation on the problem is found in 
Ref. p^ . where the instability of vortices in a binary 
mixture of BEC is studied in the weak interaction limit. 

As for the dynamical instability, another treatment is 
known, which we refer to as the RK method given by 
Rossignoli and Kowalski In this method the quan- 
tum Hamiltonian of the quadratic form of creation and 
annihilation operators is considered, and the complex 
modes appear as a result of diagonalizing it with un- 
usual operators which arc neither bosonic nor fermionic 
ones. In the previous work [21], our group has analyti- 
cally derived the condition for the existence of complex 
modes for the case where the condensate has a highly 
quantized vortex, by using the RK method. There a 
small coupling expansion is adopted, and the two-mode 
approximation on the Hamiltonian is assumed and is es- 
sential. The three- mode analysis is also performed, and 
it is confirmed that the condition for the existence of 
complex modes is not modified. However, it is not clar- 
ified why the two-mode approximation is crucial for the 
appearance of complex modes. 

In this paper, we study analytically the BdG equations 
whose two-component eigenf unctions are not only of 
positive-norm but also of negative-norm or of zero-norm. 
It is shown that the degeneracy between a positive-norm 
mode and a negative-norm one is necessary for the emer- 
gence of the complex eigenvalues. To do this, we ex- 
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pand the BdG equations in powers of a shift about a 
value of the couphng constant at which all the eigenval- 
ues are real. The analytic study in Ref. [l^ has already 
indicated the importance of the degeneracy, which was 
referred to as the frequency resonance of positive- and 
negative-energy excitations there. It is emphasized that 
the present analysis is quite simple and general, not re- 
stricted to particular systems. Furthermore, our study is 
based on the expansion of the BdG equations around a 
finite (non-zero in general) coupling constant, while the 
one in Ref. [l^ is limited to the small coupling constant. 
We apply the considerations developed here to the case 
where the condensate has a highly quantized vortex, and 
derive the expression of the condition for the emergence 
of the complex modes which completely coincides with 
our previous one [2l| . Through this analysis, the validity 
and implication of the two-mode approximation in our 
previous work is confirmed. 

This paper is organized as follows: In Section II, the 
BdG equations are introduced and their properties are 
shown. In Section III, we expand both the GP and BdG 
equations in powers of a shift of the coupling constant, 
and show that the degeneracy between the two modes 
with different norms causes the complex eigenvalues. In 
Section IV, we apply the analysis to the case where the 
condensate has a highly quantized vortex. Section V is 
devoted to the summary. 



II. BOGOLIUBOV-DE GENNES EQUATIONS 
AND THEIR PROPERTIES 



We start with the following TDGP equation to describe 
the dynamics of the condensate, 

i^^^{x,t) = [K + V{x) + g\'<SJix,t)\^]'^ix,t) , (1) 

where K = — ^JT^^ with the mass of a neutral atom 
M, and V{x) and g are the trap potential for the atoms 
and the coupling constant of self-interaction, respectively. 
Throughout this paper h is set to be unity. The quan- 
tity |^'(a;,t)p corresponds to the density of condensate 
atoms, and the total condensate number N is given by 
TV = J(Px |^'(a;,t)p. The stationary state is given as 
^{x, t) = ^(a;)e~''^* with the chemical potential fi, where 
the function ^(a;) satisfies the stationary Gross-Pitaevskii 
(GP) equation. 



[K + V{x)-^l + g\ax)\^] ax) = 0. 



(2) 



To study the excitation spectrum of the condensate, 
let us consider a small fluctuation around the station- 
ary state as 1i{x,t) = [^{x) + S^{x,t)]e~^^'*. Substi- 
tuting it into Eq. ([1]) and linearizing it with respect to 
S'i'{x,t) = M„(cc)e-'-^"* + <(a;)e*-^''*, we obtain the fol- 
lowing BdG equations. 



Here the doublet notation is introduced as 
y^{x) ^ 

T = 



where 



Vn{x) 

L M 

-M* -C 



C - K + Vix)-ti + 2g\ax)\\, 
M - ge{x). 



(4) 
(5) 



(6) 
(7) 



The eigenvalue En is not necessarily real, since T is not 
a Hermitian operator. The similar contents given in this 
section below arc presented in Ref. p^l for a general case, 
and in Refs. [13, 19, 23] for the particular systems. 

It is easily shown that the operator T has the following 
algebraic properties [H, : 



(8) 
(9) 



with the i-th Pauli matrix ai. These relations regulate 
the properties of the eigenfunctions in the following. 

Let us introduce the "inner product" for arbitrary dou- 
blets s{x) and t{x) as 

s{x),t{x)^ = jd\s\x)a^t{x), (10) 

and the squared "norm" as 

's{x),s{x)). (11) 



Note that the squared "norm" may be zero or negative 
due to the indefinite metric CT3. 

It is shown from Eq. ([S]) that the operator T is pseudo 
Hermitian, namely. 



(^six), Tt{x)^ = {^Ts{x), t{x) 



(12) 



Because of the conjugate symmetry of the "inner prod- 
uct," 



(13) 



Eq. implies that the matrix element ( s, Ts ) is 
always real. Furthermore, Eq. (fT2)) leads to 



{Er,-ElM Urn, ?/„ - 



and 



im(£;„)lly„ll' = 0, 



(14) 



(15) 



Tyn{x) = EnVnix) 



(3) 



where y„ and ym are the eigenfunctions of the BdG equa- 
tions which belong to the eigenvalues En and Em , respec- 
tively. Two eigenvectors y„ and ym are orthogonal to 
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each other for ^ Em- If En is a complex eigenvalue 
[lm{En) ^ 0], then we have |ly„|l^ — 0. Equivalently, if 
||y„|P 7^ 0, then £'„ is real. 

Next, if y is an eigenfunction belonging to a non- 
zero real eigenvalue E, one can show from Eq. ^ that 
z = aiu* is an eigenfunction belongin g to —E . One als o 

12 



where 



finds II y I 



|2;|| , as is seen in Refs 



i g to — -C/. Une also 

[ilMIillli,!!!- 



Hereafter, we denote positive- and negative-norm eigen- 
functions by y„ and Zn, respectively, which make a pair 
and are related to each other as Zn = criy*. Without 
losing the generality, the absolute values of their norms 
can be set to unity. 



Vn, Vn 



= Snr. 



(16) 
(17) 
(18) 



with 





/ Co 


Mo 
-Co 


r = 


( ^' 

\-M'* 


M' 

-a 



(28) 
(29) 



Co = i^ + F-A^W+2go|e<°f , (30) 

£' = 25'|^(")p+45oRe(e(°)e«*)-M(i), (31) 

Mo - .goC'"'' , (32) 

M' = ff'e(°)2 + 2goC(°)^W. (33) 

Note that the operator T' is pseudo-Hermitian. 

In the case where no degeneracy is present, the first- 
order eigenvalue is 



III. GENERAL CONDITION FOR EMERGENCE 
OF COMPLEX EIGENVALUES 

In this section, in order to derive the condition for the 
emergence of complex eigenvalues, we employ the pertur- 
bation theory, dividing the coupling constant g as 



(19) 



with a sufficiently small parameter e. We consider the 
situation that all the eigenvalues are real for the unper- 
turbed coupling constant go and that the complex mode 
arises due to the perturbation eg' in Eq. p^ . This ap- 
proach is quite general and is applicable to most dynam- 
ical instabilities of the BECs such as the splitting of a 
vortex, the decay of the flowing condensate in an opti- 
cal lattice. Similarly to the following analysis, it is also 
possible to develop the perturbation theory not of the 
coupling constant but of the strength of the optical lat- 
tice potential for the flowing condensate in an optical 
lattice, although we do not present it in this paper. 

The quantities appearing in the GP and BdG equa- 
tions, ([T]) and ([3]), are expanded in terms of e as 

ax) = ^^"\x)+e^^'\x) + 0{e^), (20) 

y{x) = y(")(a;)+£y(i)(a;) + 0(e2), (21) 

^ = /i(")+£/i(i)+0(e2), (22) 

E = +0(e2). (23) 

We obtain at the orders of e'~' and e^, 

CoC^"^ -Mo^^°^* =0, (24) 

Co^^^^ - Mo^'^^^* + C'C^°^ ~ M'C^°> = , (25) 

(To yW(a=) = 0, (26) 

(To - y^'Hx) + (t' - y(°)(a:) = , (27) 



(34) 



Since T' is pseudo-Hermitian, i?^^^ is always real [see the 
discussion below Eq. (fT2|) ]. Hence, if all the eigenvalues 
are real and there is no degeneracy for the unperturbed 
coupling constant go, then no complex eigenvalue exists 
for 5 = go + eg'. 

Next, let us turn to degenerate cases. Here we assume 
for simplicity that only a single pair of modes is degen- 
erate. 

First, we consider the case in which two positive-norm 

modes and yj*^-* share the same eigenvalue i?^"-* which 
is real. Using the degenerate perturbation theory, we 
obtain the following secular equation which determines 



(35) 

This implies that E^-^^^ is the eigenvalue of the Hermitian 

matrix, 



yf\T'y^' 



T'yl 



(0) 



T'y'i 



(0) 



(36) 



hence E^^^^ is real. Likewise, it is also shown that no 
complex eigenvalue arises from the degeneracy between 
two negative-norm modes, zf''' and z^"^ . 

Suppose that the two modes of one positive- and one 
negative-norms, yf''' and z^\ are degenerate for a real 
eigenvalue . We then obtain the secular equation, 



{y['\T'y^; 
{4\T'y[ 



(0) ■ 



T'zi°^ 



T'z 



(0) 



= 0, 
(37) 
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which impHes that i?^^^ is the eigenvalue of the non- 
Hcrmitian matrix 



(0) 



'-2 1 




(38) 



In this case E'--^^ can be complex. Actually, i?*^^^ becomes 
complex if the following condition is satisfied, 



^2 1 ^ -^2 



< 2 



(39) 



where Eqs. ([T^ and are employed. 

We conclude that the degeneracy between modes of 
positive- and negative-norms, y and z, is essential for 
the emergence of the complex eigenvalue. This can 
be extended to the cases where more than two modes 
are degenerate, because the degeneracy between y's and 
2;'s gives rise to a non-Hermitian matrix such as ([38]). 
Our conclusion is consistent with the suggestions in 
Refs. [l3|,[l3]. Note that for the emergence of the complex 
eigenvalues the degeneracy is necessary but not sufficient 
and that the inequality in Eq. (|39p must be satisfied in 
addition. 



IV. APPLICATION TO THE CASE OF A 
HIGHLY QUANTIZED VORTEX AND VALIDITY 
OF TWO-MODE APPROXIMATION 

In this section, as an application of the formulation 
developed in the preceding section, we derive the expres- 
sion of the condition for the emergence of the complex 
eigenvalues in a trapped BEG with a highly quantized 
vortex. The application elucidates why the two-mode 
approximation is valid in our previous work where the 
RK method is used [U. 

We consider a harmonic trap potential with a cylindri- 
cal symmetry: 



F(r,z) = iM(c.ir2+a;,V), 



(40) 



where r = ^ + y^- Let us assume that the vortex is 
created along the z-axis. Then the order parameter ^[x) 
can be written as 



(41) 



with the winding number k which is an integer. Without 
losing the generality, k is assumed to be a positive integer. 
We consider the case of go = 0, which corresponds to the 
weak coupling limit. The quantities £o, M.^, and M! 
are written as 



C 
Mo 
M' 



K + V 
,(1) - 



,(0) 



0, 



(42) 
(43) 
(44) 
(45) 



The eigenf unctions of the unperturbed BdG equations 
are written as 



AO) 



Un,l,rn{r,6,z) 









(46) 



(47) 



which belong to the eigenvalues en,i,m and — eri,;,m, re- 
spectively. Here, n (= 0, 1, 2, . . . ) is the principal quan- 
tum number, I (— 0, ±1, ±2, . . . ) is the magnetic quan- 
tum number, and m (= 0, 1, 2, . . . ) is the quantum num- 
ber along the z axis. The function Un^i^rniTT^^^z) is ex- 
pressed as 

Un^iA^, 0, z) = yTe»('+«)%„ ; „(r, z) , (48) 

with the solution Un,i,7n{f, z) of the following eigenequa- 
tion: 

'2M I 9^ + :p2 + ^ 



+V{r,z)-fi 



(0) 



Un,l,mir,z) 



(49) 



The explicit form of Un^i.miT, z) is given by the Laguerre 
polynomials L^{x) and the Hermite polynomials Hm{x) 
[24| as 

u„,i,™(r, z) = C„,z,„e-5"i'-'(a^r)l'+«lL^'=l(air2) 

X e-3"^'iJ„(a^z), (50) 

with 



7rV2 Y 2™ m! (|Z + k| +n)! 
while the eigenvalue is given as 



(51) 



en,;,m ^ uj± {2n + \l + k\ + 1) + ujz m + 



(52) 



The normalization constant (15111 is determined so that 



d xU*i^{r,0,z)Un',l',m'ir,0,z) =6„n'6w6r, 



(53) 

and||yS,„|P = -||.(°Ull' = l- 

The unperturbed order parameter and the chemical 

potential are obtained as 

(54) 
(55) 

Substituting Eqs. dUD-din]) into Eq. we obtain 

(2k)! , . 

/id) = d'-—{^ , (56) 



e(°)(r,0,z) = ViVC/o,o,o(r,0,z), 



(k!)2 ' 
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with 



Now, let us evaluate the matrix elements of T' 



D = 



(57) 



\ n,i,m' n' ,1' ,m' J 

( (0) (0) 



X 9'N f, , 2 

dl.-l' — / drdz q oUn,l,mUn' .-l,m' ■ 

In . 



(58) 
(59) 
(60) 



r 



We notice that the matrix elements in Eqs. ([58)) and 
vanish unless I = l\ while that in Eq. (|60|) vanishes unless 
I = —I'. According to the general formulation of the de- 
generate perturbation theory, the whole secular determi- 
nant is a product of secular determinants of submatrices 
corresponding to the subspaces spanned by the degener- 
ate eigenf unctions. We therefore inquire the degeneracy 
in Eqs. (HSD-dSOl)- Concerning Eqs. dHSD and dH]), the 
degeneracy condition en,i,m = ^n',i'.m' implies I — I' , and 



then n = n' and m = m! . On the other hand, concerning 
Eq. ([60)) . the condition is given by en,i.m = —^n',i',m', and 
is satisfied when I = —I' and n = \ 



< K. Consequently, only the pairs of j/q /n and 



m — m' = and 

'\,lo ^^'^ ^0,-1,0 

are degenerate and give the non-trivial secular determi- 
nants. Thus, the whole secular determinant is a product 
of the following three kinds of factors: 



J 



(i) 
(ii) 

(iii) 



( 



.(0) 



„ T'z 



^(0) J./ (0) 

^o,-i,o' ^ yo,i,o 



.(0) 



T'z, 



(0) 



^o,-;.0' ^ ^0,-1.0 1 + ^ 



It is clear that the only term which can cause com- 
plex is (iii), which is the form of the determinant of a 
2x2 matrix, as discussed in the preceding section. This 
analytic result justifies the two-mode analysis assumed 
in our previous work [2l[, and clarifies which two modes 
are relevant. 

Now the condition for the emergence of the complex 
mode is given by 



JO) 



'0.-i,0 



< 2 



ryS,)\ ■ (61) 



Calculating the inner products, we finally obtain the ex- 



pression of the condition as 



(2k + 0! 



(2k- 0! 



2'k!(k + /)! 2-'k!(k-0! 



< 



(2k) 



(k!)2 
(2k)! 



K\y^{K + iy.{K-l)l 



(62) 



We can see from the inequality that the complex eigen- 
values appear when \l\ — 2 with k = 2, \l\ =2,3 with 
K = 3, or \l\ = 2,3 with k = 4. In the case of k = 0, 1, 
the inequality is never satisfied and no complex eigen- 
value appears. This result is consistent with that of the 
numerical studies of the BdG equations [l^, Ell i ^^'^ 
produces our previous one j21| . 

In our previous analysis [21j , we did not treat the BdG 
equations but the RK Hamiltonian with the two-mode 
approximation, although its justification was not clear 
then. In this section, we have confirmed for the system 
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of a highly quantized vortex system that the two-mode 
analysis in our previous treatment is justified, as the cru- 
cial factor of the whole secular determinant is the secular 
determinant of a 2 x 2 matrix. Furthermore, it is now 
clear which two modes are relevant. 



V. SUMMARY 

We have investigated the condition for the appearance 
of the dynamical instability of the Bose-condensed sys- 
tem, whose sign is interpreted to be the emergence of 
the complex eigenvalues in the BdG equations, using the 
perturbation theory for both the GP and BdG equations 
with respect to the coupling constant. We note in com- 
parison with Ref. [l^ that our method is quite simple and 
general without being confined to particular systems, and 
that the perturbative expansion is made around a finite 
(non-zero in general) coupling constant. 

Our important conclusion is that the emergence of the 
complex eigenvalues in the BdG equations is attributed 
to the degeneracy between a positive-norm eigenmode 
and a negative-norm one. The expression of the condi- 
tion for the existence of the complex eigenvalues follows 
from the secular equation of the 2x2 matrix. We have 
considered the situation that all the eigenvalues are real 
for the unperturbed coupling constant 170, and that the 
perturbation in the coupling constant gives rise to the 
complex modes. 

Inspired by the analysis on the flowing condensate in 
an optical lattice with sufficiently small lattice height 
p^ . suggesting that the complex eigenvalues are caused 
by the degeneracy, Taylor and Zaremba solved the BdG 
equations in the weak potential limit and confirmed the 
suggestion [2^. Similarly as in the case of the fiowing 
condensate in an optical lattice, the perturbation with 
respect to the potential strength can be considered. Fol- 



lowing the similar discussion as we developed in this pa- 
per, we develop the perturbation theory not of the cou- 
pling constant but of the potential strength. Then it will 
turn out that the appropriate degeneracy at the 0-th or- 
der of the BdG equations is necessary for the emergence 
of the complex eigenvalues. In this way, our method pre- 
sented in this paper includes the method of Taylor and 
Zaremba. 

Furthermore, we have applied the analysis to the case 
where the condensate has a highly quantized vortex, and 
have rederived the expression of the condition for emer- 
gence of the complex eigenvalues. We have confirmed 
why the two-mode approximation applied in our previ- 
ous work [2l| is valid, deriving the form of the secular 
equation, which is a product of 2 x 2 matrix or less. 
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